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Abstract 

In earlier work it was shown that a weak modification of general rel- 
ativity, in the linearized approach, renders a spherically symmetric and 
stationary model of the Universe. This was due to the presence of a third 
mode of polarization in the linearized gravity in which a curvature en- 
ergy" term is present. Such term was identified as the Dark Energy of the 
Universe. 

In this letter a more reahstic model is discussed. A different cosmo- 
logical solution to the Einstein- Vlasov System is analysed. This solution 
shows reasonable results which are within the standard bounds predicted 
by the cosmological observations. 

PACS numbers: 04.50.Kd, 98.80.Jk. 

Over the past century, a standard cosmological model has emerged. With 
relatively few parameters, the model describes the evolution of the Universe 
and astronomical observations on scales ranging from a few to thousands of 
Megaparsecs. In this model the Universe is spatially fiat, homogeneous and 
isotropic on large scales, composed of radiation, ordinary matter (electrons, 
protons, neutrons and neutrinos), non-baryon Cold Dark Matter, and Dark 
Energy. The galaxies and the large-scale structures grew gravitationally from 
tiny, nearly scale-invariant adiabatic Gaussian fiuctuations [U [5] . The WMAP 
data offers a demanding quantitative test of this model [3]. In other words, the 
Universe is seen like a dynamic and thermodynamic system. The test masses 
(i.e. the "particles") are the galaxies which are stellar systems with a number 
of the order of 10^ — 10^^ stars. The galaxies are located in clusters and super 
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clusters, and observations show that, on cosmological scales, their distribution is 
uniform. These assumptions can be summarized in the Cosmological Principle: 
the Universe is homogeneous everywhere and isotropic around every point [U [5] . 
The Cosmological Principle simplifies the analysis of the large scale structure. 
It implies that the proper distances between any two galaxies is given by a 
universal scale factor which is the same for any couple of galaxies [H [5] . 

On the other hand, although the standard model of the Universe achieved 
great success, it also displayed many shortcomings and flaws. As a matter of 
fact, the accelerated expansion of the Universe, which is observed today, implies 
that cosmological dynamics is dominated by the so called Dark Energy, which 
gives a large negative pressure. This is the standard picture, in which this new 
ingredient is considered as a source on the right-hand side of the field equations. 
It should be some form of un-clustered, non-zero vacuum energy which, together 
with the clustered Dark Matter, drives the global dynamics. This is the so called 
"concordance model" (ACDM) which gives, in agreement with the CMBR, LSS 
and SNela data, a good picture of the observed Universe today, but presents 
several shortcomings such as the well known "coincidence" and "cosmological 
constant" problems [4]. An alternative approach is changing the left-hand side 
of the field equations, to see if the observed cosmic dynamics can be achieved 
by extending General Relativity [5]. In this different context, it is not required 
to find candidates for Dark Energy and Dark Matter, that, till now, have not 
been found. Only the "observed" ingredients, which are curvature and baryon 
matter, have to be taken into account. Considering this point of view, one can 
think that gravity is different at various scales and there is room for alternative 
theories. In principle, the most popular Dark Energy and Dark Matter models 
can be achieved considering f{R) Theories of Gravity, where R is the spacetime 
curvature, see the review [5], the recent results [B]-[in] and references within. In 
this picture, the sensitive detectors for gravitational waves (GWs), like bars and 
interferometers, whose data analysis recently started could, in principle, be 
important. In fact, a consistent GW astronomy will be the definitive test for 
General Relativity or, alternatively, a strong endorsement for Extended Theories 
of Gravity [H]. 

In earlier work [T3] it was shown that a weak modification of General Relativ- 
ity, in the linearized approach, renders a spherically symmetric and stationary 
model of the Universe. This was due to the presence of a third mode of polar- 
ization in the linearized gravity in which a ^'curvature energy^^ term is present. 
Such term was identified as the Dark Energy of the Universe. 

In this letter a more realistic model is discussed. A different cosmological 
solution to the Einstein- Vlasov System is analysed. This solution shows reason- 
able results which are within the standard bounds predicted by the cosmological 
observations. 

Let us consider the action 




(1) 



This is a particular choice in /(i?) Theories of Gravity [5]-|10j. 
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In the limit e — >■ and /o = 1 the action ([T|) recovers the canonical form of 
the Einstein-Hilbert action of General Relativity [TJ [2] , i.e. 



S = / d'^x^R + (2) 



Criticisms on f{R) Theories of Gravity arise from the fact that lots of such 
theories can be excluded by requirements of Cosmology and Solar System tests 
[14]. However, in the case of the action ([1]), the discrepancy with respect to 
standard General Relativity is very weak, because e is a very small real pa- 
rameter. Thus, the mentioned constraints could, in principle, be satisfied. In 
particular the authors of [13] found 

0<e< 7.2*10-1^. (3) 

Fundamental constrains can be renormalizcd in order to obtain /o = 1 in the 
action ([T]) 

Wc analyse interactions at cosmological scales, therefore the linearized the- 
ory in vacuum, i.e. with Cm = 0, can be considered. It gives a better approxi- 
mation than the Newtonian theory [I] [2] [13]. Thus, let us start from the pure 
curvature action 

S = j d^x^gUR^+'. (4) 

The theory arising from such an action has been linearized in |13| . A short 
review of the linearized approach is needed for a better understanding of the 
theoretical framework. 

By varying the action ([4]) with respect to g^j^ , the field equations arc obtained 
(through this letter the convention G = 1, c = 1 and h ^ \ will be used and 
Greek indices run from to 3)) [13] 

G,. - ^^^l^J^^^ -\9^..ehR'+' + [{l+e)hR%,,,-g^,U[{l+e)hR']}. (5) 

The trace of these field equations is 

□ (1 + e)hR' = ^l^foR'+^. (6) 
By making the identifications |13j 

^^il+e)foR- and ^ ^ (7) 

a Klein - Gordon equation for the effective $ scalar field is obtained. It can 
be written as 

,s, 

Let us consider a little perturbation of the background, which is assumed 
given by a near Minkowskian background, i.e. a Minkowskian background plus 
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$ = $0 (the Ricci scalar is assumed constant in the background) [12]. We also 
assume $0 to be a minimum for the effective potential V: 



V ~ ia(5$2 ^^c^ EH^, (9) 

and the constant E has mass-energy dimension. 
Let us write 

(10) 

$ = $0 + (5$. 

To first order in ft,^^ and (5$, calling Rfj,^pa , and R the linearized 
quantity which correspond to R^^pa , Rfj.^ and i?, the linearized field equations 
are [13] 

(11) 

where 

Notice that even if the minimum of the effective potential V is the real 
number 0, V is not equal to $, therefore we are not dividing by in Eq. (|T2]) . 
In other words, as ^ 7^ $, the condition V = Q does not imply $0 = 0. 

Then, from the second of Eqs. (jH]), one can define the mass-energy like [13] 



Thus, as the mass-energy is generated by variation of the Ricci scalar, we can 
say that, in a certain sense, it is generated by variation of spacetime curvature. 
This is exactly the curvature energy''^ term |13j . 

Wc recall that, in the present case, the theory is viable as the modifica- 
tion of General Relativity is very weak and in agreement with requirements of 
Cosmology and Solar System tests |14| . 

In |13| it has been shown that, by using gauge transformations, the lineariza- 
tion process gives 

V(i, = A+{t - z)eW + AX (t - z)eS^^J + h,{t, z)ijp,. (14) 

The term A+(t — z)e\X' + A^{t — z)e|i^^ describes the two standard polar- 
izations of gravitational waves which arise from General Relativity, while the 
term hg(t, z)r]fj^i, is the massive field arising from the high order theory. In 
other words, the function R^ of the Ricci scalar generates a third polarization 
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state for GWs which is not present in standard General Relativity. This third 
polarization has a "curvature" energy E, see [13] for details. 

In the simple model of Universe in |13| the dynamic of the matter is described 
by the Vlasov equation [T31 [TS]. The gravitational forces between the particles, 
i.e., the galaxies, are supposed to be mediated by the third mode of Eq. ^(T^ 
after making the assumption that at cosmological scales such a mode becomes 
dominant (i.e. A~^,A~ <C h^) [IS]- In this way the "curvature" mass-energy E 
can be identified as the Dark Energy of the Universe ~ lO^^^g/cm^ [5]. 

The model in |13j realized a spherically symmetric and stationary model of 
the Universe. The present work analyses a different cosmological solution to 
the Einstein- Vlasov System. The new model shows reasonable results which are 
within the standard bounds predicted by the cosmological observations. 

In the hypothesis A^,A~ <^ h^, Eq. (|14p can be rewritten as 

hf^i,{t, z) = he{t, z)rif_^^, (15) 
and the line element of the model is the conformally flat one 

ds^ = [1 + h,{t, z)]{dt^ - dz^ - dx^ - dy^). (16) 

By defining 

a' = l + K(t.z), (17) 

the line element ([T6| results similar to the cosmological Fricdmann- Robertson- 
Walker (FRW) line element of the standard homogeneous, isotropic and flat 
Universe [U E] 

ds^ = [a^{t, z)]{+dt^ - dz^ - dx^ - dy^). (18) 

Strictly speaking, this metric does not describe an homogeneous and isotropic 
universe, as is also a function of z. Thus, we have to further assume dji^ = 0, 
which removes the z-dependence [TB] . 

To satisfy the condition demanding that the particles make up an ensemble 
with no collisions in the spacetime, the particle density must be a solution of 
the Vlasov equation 

ft/ + ^a.»/-r;i,^9,./ = o. (19) 

Here F^^ represent the usual connections, / is the particle density and is 
determined by p°'{a = 1, 2, 3) according to the relation 

5m-P>"^ = -1 (20) 

[131 115] , which expresses the condition that the four momentum p'' lies on the 
mass shell of the metric. 

A clarification is needed. We have to distinguish between the standard ki- 
netic theory which is based upon Boltzmann equation [17j and the Einstein- Vlasov 
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system which we are discussing in the present work. In a kinetic description the 
time evolution of the system is determined by the interactions between the par- 
ticles, which depend on the physical situation [TS]. For instance, the driving 
mechanism for the time evolution of a neutral gas is the collision between parti- 
cles (the Boltzmann equation) [T5] . For a plasma the interaction is through the 
electromagnetic field produced by the charges (the Vlasov-Maxwell system), and 
in astrophysics the interaction between collisionlcss particles is gravitational (the 
Vlasov-Poisson system and the Einstein- Vlasov system) [TB]. Therefore, an im- 
portant difference between the standard kinetic theory and the Einstein- Vlasov 
system which we are discussing in this letter is that in the first case there are col- 
lisions between particles while in the Einstein- Vlasov system particles interact 
gravitationally only and they are collisionless. 

We recall that, in general, the Vlasov- Poisson system is [131 [15] 

dtf + V ■ Vxf - \/xU ■ \/.af = 



AC/ = 47rp (21) 

p{t,x) = J dvf{t,x,v), 

where t denotes the time and x and v the position and the velocity of the 
galaxies. The function U = U{t,x) is the average Newtonian potential gener- 
ated by the galaxies. This system represents the non-relativistic kinetic model 
for an ensemble of particles with no collisions, which interacts through the grav- 
itational forces that they generate collectively [131 [T5|. Thus, one can use such 
a system to describe the motion of galaxies within the Universe, thought of as 
pointlike particles, when the relativistic effects are negligible [131 [15]. In this ap- 
proach, the function f{t, x, v) in the Vlasov- Poisson system pip is non- negative 
and gives the density on phase space of the galaxies within the Universe. 

By using the classical transformation from conformal time to synchronous 
time [1] 

dt -> '^^ (22) 

the line element p6l) . in spherical coordinates, becomes 

ds^ = dt^ -[l + K {t)] (dr^ + {dO^ + sin^ Odip^ ) ) . (23) 
The metric tensor has the form 




= I r (24) 

where 7m„ l + h^{t). 



6 



Following [2], we define Xmn = ■§["fmn. The Einstein field equations in the 
synchronous frame are: 

R°o = - Ixlxt = (T° - It) (25) 

K = l{xU-xt,a)=T^ (26) 

R'a = -P'a ^§iiV^xl) = - \5'aT, (27) 
where is the Ricci tensor in 3 dimensions [2]- 

On the other hand, the Einstein field equations in the Einstein- Vlasov system 
are [Ml US] 

G,,, = j fit, x,p)p^pJip'' + m'')d^p, (28) 

where m is the mass of a particle (galaxy). 

We can split the function f{t,x,p) into a couple of equations for f+{t,x,p) 
and f-(t,x,p) which are constructed by reducing f(t,x,p) respectively on the 
"upper" half and on the "lower" half of the mass shell. Eq. (|T9)) becomes 



dtf± = -—{ 7">n -^PnPr j^)f±- (29) 

p'i \ dx,n 2 dXm Op,n J 

Eq. pop can be interpreted in Hamiltonian terms: 

pldtf±^{HJ±}, (30) 
where the Hamiltonian function is 

H ^ \l"'-praPn. (31) 

One can calculate the components of energy-momentum tensor T^i, in the 
approximation which considers galaxies like massless particles (m = in Eq. 

I) 



' IJ±±i^Ji±^^^S, (32) 



T — 



{^l + K{t)fr'^ sin 61 y yiTM^V r2 r2sin6l" 

^ ^ .JT+hAt)^±fiI^=p^p^d'p (33) 



sin t 



{^l + K{t)fr^ sm0 J " J p'^+pI I pI 

1 

{^\ + K{t)fr'^ sin^ 



To™ = -^=_^— — ^ / (/+ - f^)p,nd'p. (34) 
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The Einstein field equations (pSj). ([^ and (P7)) give two independent dy- 
namic equations which can be written down in terms of the scale factor a = 

a^^-l + ^ f (/+(,s) + f-(s))-d^s (35) 
6a J a 

■ ^ 2- + \ ({U{s) + f^{s))d^s, (36) 
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a a 



where 



2 2 

2 I f^2 I -Pa /'n>7\ 



By introducing the dimensionless variables r and t we put 

a = aor 

i = lot. 

• dr 

-~ It 



(38) 



where oq is the present value of the scale factor of the Universe. 
Eq. ([55]) becomes 

f = -1 + ^ 

(39) 

210 = 1- 

The solution of the system ([39|) is 



z:(i)-V^-(^-vJ^)' (40) 

if j > 1. Notice that expression PO)) becomes imaginary for certain values of 
the parameter t [17j . Such a parameter represents the cosmic time normalized 
by the present value of the scale factor ao , see Eq. ([38)) . By choosing the origin 
of the cosmic time at the present era of the cosmological evolution, i.e. to = 0, 
expression (|40p is always real if j > 1. 

Returning to the {t, a) variables we get: 



.2_ ^ 1^2, 



The today's observed Hubble radius and the today's observed density of the 
Universe arc respectively ao lO^^cm and po sa 10~^^cm~^. Therefore 

By inserting these data in Eq. (|4ip we obtain a singularity at a time 
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ts ~ —10^° seconds and a value for the today's theoretical Hubble constant 
Ha = 10-29cTO-i. 

As it follows from the above analysis, even under the assumption to neglect 
the baryon mass of the galaxies the results look reasonable. They are of the 
same order of magnitude of the standard cosmological model [3]. 

Conclusion remarks 

In earlier work it was shown that a weak modification of general relativity, in 
the linearized approach, renders a spherically symmetric and stationary model 
of the Universe. This was due to the presence of a third mode of polarization 
in the linearized gravity in which a curvature energy^^ term is present. Such 
term was identified as the Dark Energy of the Universe. 

In this letter a more realistic model has been discussed. A different cosmo- 
logical solution to the Einstein- Vlasov System has been analysed. This solution 
shows reasonable results which are within the standard bounds predicted by the 
cosmological observations. 

We hope in further analyses which could insert the baryon mass too and 
realize a better fine-tuning the model with the cosmological observations. 
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